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Abstract: We present a conceptually new framework for describing jet evolution in 
the dense medium produced in ultra-relativistic nucleus-nucleus collisions using pertur- 
i bative QCD and its implementation into the Monte Carlo event generator Jewel. The 

rescattering of hard partons in the medium is modelled by infrared continued pQCD ma- 
J> trix elements supplemented with parton showers. The latter approximate higher order 

real-emission matrix elements and thus generate medium-induced gluon emissions. The 
interplay between different emissions is governed by their formation times. The destruc- 
tive interference between subsequent scattering processes, the non-Abelian version of the 
Landau-Pomeranchuk-Migdal effect, is also taken into account. In this way the complete 

£S) radiation pattern is consistently treated in a uniform way. Results obtained within this 

minimal and theoretically well constrained framework are compared with a variety of ex- 
perimental data susceptible to jet-quenching effects at both RHIC and the LHC. Overall, 

^ a good agreement between data and simulation is found. This new framework also allows 

to identify and quantify the dominant uncertainties in the simulation, and we show some 
relevant examples for this. 
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1 Introduction 

High momentum transfer processes occur abundantly in hadronic interactions at collider 
energies. This is not only true for proton-proton collisions at the LHC and other similar 
collider experiments before it, but also for nucleus- nucleus collisions at RHIC and at the 
LHC. Such processes are well understood and calculable in perturbative QCD and they can 
be simulated with modern Monte Carlo event generators [1-3] using QCD matrix elements 
supplemented with parton showers [4] . The latter resum the dominant collinear logarithms 
associated to the emission of additional partons to leading (and partly sub-leading) loga- 
rithmic accuracy, with many improvements and reformulations in the past decades [5-13]. 
In the present work, we discuss how such techniques can be used to formulate a general 
framework for the parton shower evolution in the presence of dense QCD matter produced 
in nucleus- nucleus collisions. 

Simple considerations, based on scale arguments and the uncertainty relation, indicate 
that final state parton showers evolve for up to several fm/c, which is comparable to the 
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transverse size of the dense QCD matter produced in the spatially extended overlap region 
of ultra-relativistic nucleus-nucleus collisions. Therefore it can be expected that the QCD 
parton showers introduced and studied in the vacuum will be subjected to medium modifi- 
cations in nucleus-nucleus collisions. Data from RHIC and the LHC strongly support this 
view. In particular, in comparison to standard perturbative benchmarks from pp-collisions, 
single inclusive hadron spectra at high transverse momentum are strongly suppressed, by 
up to a factor of 5 at RHIC [14, 15], and by up to a factor of 7 at the LHC [16-18]. 
Over the last two years, data on reconstructed jets in heavy-ion collisions have provided 
further, more differential information about such medium modifications of QCD radiation. 
In particular, measurements of the di-jet asymmetry [19, 20] and the inclusive jet suppres- 
sion [21, 22] indicate that the interaction of parton showers with the surrounding QCD 
medium can dissipate a significant fraction of its total transverse energy to large angles, 
outside typical jet cones . On the other hand, the fragmentation functions inside the jet 
cone appear to experience only moderate medium modifications [25, 26]. 

The modification of leading fragments has been anticipated more than two decades 
ago [27-29] and triggered the development of a large number of jet quenching models 
in subsequent years. For recent reviews on such models and their comparison to data 
from RHIC and the LHC, see refs. [30-33]. In the following, we limit our discussion to 
highlighting the main similarities and the main differences between our proposal and some 
of these previous works. We start by recalling that there are essentially three qualitatively 
different analytically formulated approaches to medium-induced parton energy loss, all of 
which based on strong assumptions about kinematics and/or medium properties. 

Firstly, starting with refs. [27, 28] there are analytical studies of collisional parton 
energy loss, in which highly energetic partons lose energy and randomise momentum in 
consecutive partonic 2—^2 scatterings with partonic components of the medium. These 
studies do not assume particular kinematical constraints. They are strongly limited, how- 
ever, in that they do not include medium-induced parton emission processes. 

A complementary mechanism, based on analytical studies of radiative parton energy 
loss through modifications of the emission pattern of secondary partons, was proposed in 
refs. [34-39]. For technical reasons, the emissions in this approach have been limited to 
close-to-eikonal kinematics, neglecting recoil effects. Also, typically the emission of only 
one gluon from a single highly energetic parton is considered. Adjusted with a number of 
phenomenologically motivated requirements, the simple picture introduced in the publica- 
tions above has been employed in model studies in the past years [40-43] . Recent analytical 
studies improve on some approximations by considering gluon emission from eikonal anten- 
nae [44, 45], and by providing a differential formulation of the colour flow in the multiple 
scattering process [46]. 

In addition, radiative parton energy loss calculations illustrated explicitly that by 
scattering on several spatially separated partons, the distribution of emitted gluons is 

1 Reconstructed jets in heavy ion collisions are by now typically identified through the anti-fcx jet 
algorithm [23]. The jet catchment areas that the anti-fci. algorithm associates to jets are approximately 
circular[24], and we shall refer to them for simplicity as 'cones' in the following. 
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modified by destructive interference effects. In the eikonal limit, these interference patterns 
are understood in terms of the non-Abelian Landau-Pomeranchuk-Migdal (LPM) effect. 
This effect can be implemented in a probabilistic formulation based on formation time 
constraints, and lead to the first simulations taking this effect consistently into account [47]. 

By now, however, it is widely accepted that analytical studies of radiative parton 
energy loss provide only very limited constraints on phenomenological models implemented 
in simulation codes [48], mainly because they are restricted to a small part of the relevant 
kinematics and because they only insufficiently constrain how medium-effects modify the 
vacuum structure of the parton shower. There have been strong efforts in recent years 
to overcome these deficiencies in Monte Carlo formulations of medium-modified parton 
showers [49-56] . Many of these models use parton formation times to specify destructive 
interference in parton emissions and the spatio-temporal embedding of the parton shower in 
the medium. This is closely related to the procedure that we propose in the present work. 
On the other hand, the modelling of medium-induced gluon emissions typically makes 
use of analytical results obtained in the eikonal limit. Thus, although the Monte Carlo 
models can improve on aspects like implementing exact energy-momentum conservation, it 
should be stressed that they still suffer from the numerically large systematic uncertainties 
and conceptual limitations related to using the analytical results outside their strict region 
of validity. This is why we believe it is important not only to formulate models with 
tractable assumptions but also to supplement them with a well-defined way to quantify 
related uncertainties. 

Finally, analytical studies of parton energy loss in strongly coupled non-Abelian plas- 
mas with known gravity duals should be mentioned, for a recent review cf. [57]. This 
approach is qualitatively different from the one proposed here and we do not discuss it 
further in the present work. 

In this paper, a conceptually new approach to jet quenching is formulated, which con- 
sistently includes both collisional and radiative parton energy loss mechanisms. It is based 
on known properties of perturbative QCD in non-eikonal kinematics and their continuation 
to the infrared region. This formulation consistently treats the scattering of an energetic 
parton off a medium constituent in perturbative QCD using standard techniques. The 
resulting interaction is described at leading order through standard 2 — > 2 QCD matrix 
elements, which are valid in the full perturbative regime and can easily be extended into the 
non-perturbative region with some simple, predictive parametrisations. Further, following 
the standard description of all hard scattering events at colliders, large logarithms associ- 
ated to collinear singularities occurring in matrix elements for subsequent real emissions 
are resummed by the parton shower and thus dealt with in a probabilistic language. This 
picture provides a well understood and systematically improvable approximation to extra 
gluon emissions. 

In this way both collisional and radiative energy loss are generated consistently as part 
of a more general description of scattering processes in perturbative QCD. This approach 
includes in a natural way also radiation off the scattering centre and generates configura- 
tions with any number of additional partons. So-called flavour conversion processes that 
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transfer essentially the entire energy of a parton to a different parton are also automatically 
included and the framework can easily accommodate heavy quark flavours. 

While this approach offers a general, flexible and largely well controlled framework for 
describing jet evolution in the presence of a dense, strongly interacting medium, it is not 
free of assumptions: 

1. The medium as seen by the jet is a collection of partons with a certain 
distribution in phase space. For scattering processes with a very high scale 
(resolution) this is certainly a reasonable assumption. Hard processes happen on very 
short time scales, which are much shorter than any interaction among the partons 
forming the medium. For such processes, the partonic constituents of the medium 
are resolved by the jet, the scattering partner is thus a quasi-free parton. This 
essentially is the same argument that is employed, for example, in deeply virtual 
scattering and corresponding factorisation theorems. The actual assumption here 
is that this factorisation is valid for all the interactions of the jet in the 
medium. The actual phase space distribution of scattering partners is irrelevant for 
the qualitative picture underlying the formulation of jet-medium interactions, but of 
course quantitatively impact on the results. 

2. By continuing the pQCD matrix elements into the infra-red region the 
dominant effect of soft scatterings can be included. Since it is impossible to 
restrict the framework in a meaningful way to include perturbatively hard scatterings 
in the medium only, some modelling of soft scatterings is unavoidable. While clearly 
other approaches are conceivable and may be equally legitimate, a suitable infra-red 
continued version of perturbation theory could be considered as a minimal model, 
in which the exact form of the infra-red continuation and the exact choice of the 
regulator are a source of uncertainties. These uncertainties, however, can be studied 
systematically and, with sufficient data, most likely be further constrained. It should 
be noted, though, that in particular the total scattering cross section in the medium 
depends on these choices. 

3. The interplay between different sources of radiation is governed by the 
formation times of the emissions. As it is clearly unphysical to assume that 
additional emissions happen instantaneously, there will be situations where several 
radiation processes occur at the same time. In the framework presented here, each 
emission is ascribed a formation time and for two competing emissions the one with 
the shorter formation time is accepted, while the other will not be realised. However 
plausible this seems, this nevertheless goes beyond standard parton showers, which 
have no notion of space-time evolution. Furthermore, there are only parametric 
estimates of the formation time. The resulting freedom of an exact choice will yield 
different admixtures of different processes in the same situations. This, again, gives 
rise to some uncertainties, which, however, are easily quantifiable. 

4. The physical picture behind the LPM-effect derived from results in the 
eikonal limit is valid also in non-eikonal kinematics. In calculations carried 
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out in the eikonal limit it was found that there is a destructive interference between 
gluon emissions induced by subsequent scattering processes when the distance be- 
tween the scattering centres is shorter than the gluon formation time. This is the 
non-Abelian analogue of the Landau-Pomeranchuk-Migdal (LPM) effect. A proba- 
bilistic interpretation of this phenomenon based on a physical picture was derived 
and validated in the eikonal limit. The validity of this picture also in non-eikonal 
kinematics is only assumed. A straightforward generalisation of the algorithm can 
be used to include the LPM-effect in the framework presented here. 

In the remainder of this publication, the construction paradigms above will be further 
worked out. To put the different ingredients into perspective, to further develop the lan- 
guage employed, and to define our notation, however, the reader is first reminded of a few 
concepts that will be used throughout, cf. section 2. Details of the implementation relevant 
for a better understanding of our model are presented in section 3. In section 4 we con- 
front this new model with a broad range of data, starting from those vacuum observables 
in e~e + and pp collisions, which fix the parton shower and hadronisation parameters and 
thus leaves us with a minimal set of relevant parameters. By comparing to a number of 
data from RHIC and LHC, which show jet quenching effects, we highlight the versatility 
of our model. We also make a non-trivial, testable prediction there. Finally, we summarise 
with some concluding remarks in section 5. 

2 Parton shower in a medium 

In this section, the ingredients necessary to construct an in-medium parton shower following 
the paradigms outlined above will be discussed. In section 2.1 the parton shower picture 
in the vacuum is briefly summarised. The rescattering of a hard parton in the medium is 
discussed in section 2.2, and in section 2.3 the probabilistic interpretation of the LPM-effect 
in the eikonal limit is recapitulated. Details of its generalisation to the full phase space are 
given in section 3.4. 

2.1 Reminder: parton showers 

In the collinear limit, QCD real emission matrix elements and their phase space factorise 
such that the differential cross section for any process with an additional parton in the 
final state, d<r n +i, can be expressed by the cross section for the production of an n-parton 
final state modified by process-independent terms: 



Here z is the energy or light-cone momentum fraction taken by the outgoing parton, and 
the splitting kernel Pb a {z) typically is given by the (unregularised) Altarelli-Parisi splitting 
function or kernels that reduce to it in the collinear limit. The parameter t ~ k\ ~ Q 2 ~ 6 2 
is a measure for the hardness of the additional emission, to leading logarithmic accuracy this 
can be the square of the transverse momentum k± , the virtual mass Q or the emission angle 
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9. The scale /i 2 at which the strong coupling a s is evaluated is usually given by fi 2 = k 2 ^. 
While Eq. (2.1) captures all leading collinear and soft-collinear logarithms irrespective of 
the precise definition of t, the choice t = 6 22 and the identification fi 2 = k\ ensure that 
also the next-to leading logarithmic contributions are properly taken into account. 

Closer inspection of Eq. (2.1) reveals that the cross section diverges for t — > 0, as 
expected. However, emissions at low scales t of the order of a few Aqcd will typically not 
yield any resolvable parton, instead partons radiated at such scales will most likely end up 
in the same hadron as their emitter. This justifies the introduction of an infrared cut-off 
t c , usually chosen such that the transverse momentum is cut at kj_(t c ) ~ lGeV 2 . This 
cut-off of course also implies a cut-off in z thus shielding the divergent structures inside 
the splitting kernels P. 

Iterating Eq. (2.1) for any number k of additional emissions with these choices, such 

that 



thus resums, to all orders, terms of the parametric form a s [L + L ], where L denotes 
large logarithms of ratios of the scales, to/t c . It should be noted that this accuracy is 
achieved in the limit, where the number of colours is infinite, N c — > oo. Configurations 
that are of sub- leading colour are typically suppressed with l/iV 2 and come at order a^L 21 ^ 1 



The picture developed so far treats emissions as independent from each other (inter- 
ferences between subsequent emissions are sub-leading) and independent from the hard 
process, apart from the strict ordering th = to > h > *2 > • • • indicated by the G-function 
in Eq. (2.2) and apart from the exact definition of the scale th = to, which indeed de- 
pends on the process. However, it is worth noting that the choices of t ~ 9 2 and t ~ k\ 
emerge as better suited when inspecting the radiation pattern of two subsequent emissions 
since they better capture effects related to the emergence of destructive interference in 
non-angular-ordered emissions. 

The picture above can be cast into a probabilistic form, optimally suited for detailed 
simulation by introducing the Sudakov form factor, 



It can be interpreted as the probability that a parton a emits no resolvable radiation be- 
tween the starting scale th and the cut-off t c . The independence of the emissions guarantees 
that the resulting simulation, the parton shower, can be represented as a Markov chain, 
i.e. the probability for any further emission depends only on the current state of the parton 
ensemble and not on the history leading to it. For the generation of radiation off final 

2 In dipole-like showers, also the choice t — k\ appears to allow for a resummation of such next-to 
leading logarithms [58]. 
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state particles the parton shower starts at the scale of the hard core process, each further 
emissions has to be softer than the previous until the infrared cut-off is reached. 

The fact that the result of the evolution is fixed by the hard process at the higher 
and by the known initial state, the incident hadron, at the lower scale complicates the 
generation of emissions off initial state particles. It would thus be very inefficient to evolve 
from the soft, hadronic scale up to the hard scale and discard all histories that do not 
match the hard process. Instead, initial state partons are evolved backwards starting from 
the hard process with the parton density at each evolution step constrained by the parton 
distribution functions (PDFs). This modifies the Sudakov form factor such that 

Si IS \t h ,t c ,x) = ^\s a (t h ,t c ). (2.4) 

Written in such a form, with suitable Sudakov form factors, the parton shower correctly 
resums all leading and dominant sub-leading logarithms. Various methods to further im- 
prove the accuracy of the parton shower away from the soft and collinear limits of secondary 
parton radiation exist. However, since here we are mainly concerned with the properties 
of jets, it is sufficient to stress that the parton shower, due to its logarithmic accuracy, 
correctly describes the bulk of all QCD events, and in particular the properties of jets, 
which are characterised by multiple soft and collinear emissions. 



2.2 Scattering in a medium and parton showers 

To analyse how medium effects enter the picture, consider a single scattering of an energetic 
projectile off a parton in the medium at a time that is late enough such that the parton 
shower evolution of the initial hard scattering process is finished and all emerging final state 
partons are on their mass shell. This process fundamentally is not different from the hard 
initial process that produced the hard parton in the first place. Only the centre-of-mass 
energy and thus the scale of the scattering are, on average, diminished and since there is 
no natural infrared cut-off as in the case of the initial hard scatter (the criterion that the 
jets should have at least a certain p±) one has to specify how the infrared region is to be 
treated. In Jewel an infrared continuation of the matrix elements is introduced, which 
will be discussed in section 3. Then, just as in the case of the initial hard jet production 
process, leading order 2 — > 2 matrix elements supplemented with parton showers provide 
a powerful approximation to higher-order real emission matrix elements. 

The only difference with respect to the case already discussed in section 2.1 relevant 
for the implementation of such a picture is that the incident partons in such a secondary 
scatter are not constituents of an incoming nucleon. Therefore, using nucleon PDFs in the 
initial state parton shower is inappropriate and instead 'partonic PDFs' must be invoked to 
properly account for the fact that in the backwards evolution of the initial state shower the 
incident particles are known. These parton PDFs are constructed to describe the densities 
of partons in partons, by taking into account possible radiation above the resolution scale 
t c . At and below t c no radiation is possible and the parton has no associated parton 
density. The partonic PDFs are thus computed by integrating the DGLAP equations with 
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the boundary conditions 
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In this way 2 — > 2 configurations without additional radiation, that would be classified 
as elastic processes, and configurations with 3 and more final state particles, that would 
be regarded as inelastic processes or scattering with medium-induced Bremsstrahlung are 
generated at the same time, with no additional modelling bias. It should also be noted 
that matrix elements for different multiplicities of final state partons cannot be naively 
combined due to double counting issues [59], which are typically rarely addressed in the 
phenomenology of jet quenching effects. More importantly, though, these processes all 
come with the (leading log) correct relative rates (the total scattering rate is defined by 
the LO cross section as the parton shower is unitary). In principle, hard emissions could 
also be corrected to the full matrix elements using merging prescriptions such as [59, 
60], but for the purpose of this study the accuracy of the parton shower is sufficient. A 
further advantage of this description is that it is not restricted to eikonal or close-to-eikonal 
kinematics (coherence issues will be discussed in section 3.4). 

A complication arises when a rescattering occurs before the parton shower evolution 
of the previous scattering (initial or rescattering in the medium) has terminated. The 
interplay of such different sources of radiation in Jewel is governed by formation times. 
The formation time associated with an emission parametrically is of the form 



where E is the emitting parton's energy 3 . In the case of rescattering during the parton 
shower evolution, it is the emission with the shorter formation time that is formed, while the 
other emission will be discarded. This means in turn that sufficiently hard rescatterings 
reset the parton shower and restart it at a scale given by their kinematics, while soft 
rescatterings cannot initiate further radiation. 

To further illustrate these points, consider a situation where a parton in the parton 
shower (it does not matter which scattering initiated this parton shower) is going to split 

(s) 

at a scale t\ on a timescale n, while a further scattering with scale t 2 just happens 

(s) 

before that. Then a hypothetical parton shower with starting scale t 2 is initiated. If the 
hypothetical parton shower does not radiate, then the original splitting at t\ proceeds as 
foreseen. If, on the other hand, the hypothetical parton shower would initiate radiation 
at a scale t2 with corresponding formation time T2, then there are two possible outcomes: 
If T2 > 7~i the hypothetical shower is discarded and only the splitting at t\ happens. 
In the opposite case T2 < t\ the original splitting at t\ is discarded, the hypothetical 
shower becomes real and replaces the original one. The splitting at T2 takes place and 
any further radiation from the new parton shower proceeds, provided it is not disturbed 

3 This is parametrically equivalent to the form r = 2aj/fc5_, familiar from analytical calculations of 
coherent Bremsstrahlung in the eikonal limit. 
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by further rescattering. This formation time prescription is on average equivalent to the 
statement that only rescatterings which are harder than the current parton shower scale can 
resolve the virtual parton. As a consequence, individual emissions cannot be unambiguously 
associated with a certain scattering. 

2.3 Reminder: probabilistic interpretation of the LPM effect 

Analytical calculations of medium induced gluon radiation in the eikonal limit found that 
the non-Abelian analogue of the Landau-Pomeranchuk-Migdal (LPM) effect plays an im- 
portant role. This is a destructive interference between emissions initiated by subsequent 
scattering processes that occurs when the formation times of the individual emissions over- 
lap. In [47] a probabilistic formulation suitable for Monte Carlo implementation was de- 
rived, the findings relevant for this discussion will be summarised here. 

In the eikonal limit, the radiating parton has asymptotically high energy E. The 
medium constituents then appear as static scattering centres and the momentum transfer 
between the projectile and the scattering centres is purely transverse. The original BDMPS 
approach [34] and similar other approaches [35-39] operate in a kinematical regime where 
E S> oj S> kj_,qj_, where k± denotes the transverse momentum of the gluon and q± the 
momentum transfer from the scattering centre 4 . Gluon radiation is described by a Gunion- 
Bertsch cross section of the form 



For the sake of deriving and validating a probabilistic interpretation of gluon radiation 
and in particular the LPM-effect it is entirely sufficient to keep only the singular part (i.e. 
set kj_ = qj_) and ascribe all elastic rescatterings to the gluon (as only the transverse 
momentum of the gluon relative to the projectile is relevant). 

A radiated gluon decoheres from the projectile when it has accumulated a relative 
phase 



of order unity. Here, AL is the distance it travelled. From this condition the gluon 
formation time r = 2ujk\ can be read off. 

Inspection of the calculations reveals that the scatterings occurring during the gluon 
formation time act coherently to radiate the gluon. In the gluon radiation term only the 
vector sum Qj_ = ^ qj_^ of the individual momentum transfers appears. 

The probabilistic picture is that whenever an additional scattering occurs during the 
formation time, the respective momentum transfer is coherently added. There are two 
ways how this can be achieved in practice: The first one is to update the gluon phase 
accumulated up to the last scattering centre based on the gluon momentum before the 
last momentum transfer. Then the total Qj_ including the last scattering is computed and 

4 The BDMPS-Z formula for parton energy loss also holds in the kinematical region lj ^> E — lj 3> k±, q±, 
but an explicit derivation without strong energy ordering is not known. We thank P. Arnold for this 
comment. 
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the gluon k± is changed accordingly. The new formation time is determined taking the 
already existing phase into account. The other option takes the multiple scattering during 
the formation time less literally. In this case, whenever a further momentum transfer is 
added, the new and k± are computed and the gluon formation time given by the 
updated k± is calculated. If the last scattering still is within the updated formation time 
limit (counted from the first scattering) it is accepted, otherwise it is rejected as it leads 
to an inconsistent configuration. Finally, in both cases, when there are no more coherent 
momentum transfers the gluon emission is accepted with probability l/iV sca t, where iV sca t 
is the number of coherent scatterings. 

3 Details of the implementation 

3.1 Event generation 

In our model, event generation proceeds as follows: 

1. The hard matrix elements initially producing the di-jets and the corresponding initial 
state parton shower are generated by Pythia 6.4 [61] running with the virtuality- 
ordered parton shower to provide the best fit to the Jewel simulation. There are no 
multi-parton interactions and the Cteq6l1 [62] PDFs as provided by Lhapdf [63] 
are used; simulating nuclear collisions the nuclear modification of Eps09lo [64] is 
employed. 

2. Then Jewel takes over and selects the impact parameter of the event according to the 
geometrical cross section and the transverse position of the hard scattering according 
to the density of binary nucleon-nucleon interactions in the transverse plane. It then 
generates the final state parton showers including interactions in the medium. This 
includes rescatterings, emission of secondary partons, an implementation of the LPM 
effects and so on. These details will all be discussed below. 

3. Finally, the event is handed back to Pythia for hadronisation. 

Note that in the absence of a medium the procedure is the same and the Jewel parton 
shower becomes an ordinary vacuum parton shower. This option is used for the validation 
of the Jewel parton shower, for the generation of the p+p baseline, and for the fixing of 
shower and hadronisation parameters. 

In the following we will highlight aspects of phase 2 of the sketch above. 

3.2 In medium matrix elements and their infrared continuation 

The interaction of the parton shower with the partonic constituents of the medium depends 
on properties of the medium. We specify the cross sections for 2 — > 2 processes as 



\i\mex(E,T) X. 



max 
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where the PDF takes into account possible initial state radiation off the energetic projectile. 
Note that here, implicitly, we neglect a similar evolution experienced by the target, i. e. 
initial state radiation off the medium parton. This choice, although inconsistent at first 
glance, is motivated by the fact that these partons have only very small energy in the 
laboratory frame, which translates into them being virtually unable to emit any resolvable 
parton. 

The maximum momentum transfer |t| max is determined by the initial kinematics of 
the scattering. Neglecting the scattering centre's momentum |i| max = 2m s (T)[E p — m p ], 
where m s (T) stands for the (temperature dependent) scattering centre's mass and E p and 
m p are the projectile parton's energy and (virtual) mass, respectively. The boundaries 
on the ^-integral are obtained from the requirement that k\ > (Jo/4 and are given by 
^min(|i|) = Qo/(4|t|) and x max (|£|) = 1 — Qo/(4|i|). For the partonic cross section we keep 
leading terms only, but we regularise them with a Debye mass hd ~ 3T. They therefore 
read 

%(>, lil) - ft J^dfl + A)^^- - M|f| + A)j^ ■ (3-2) 



3.3 Parton shower 

The Jewel parton shower [51] is a virtuality-ordered monopole shower, similar to the 
Pythia 6 virtuality ordered shower [61], but without applying, "by hand" the a posteriori 
angular ordering constraint, when being run in the medium. Identifying the evolution 
parameter t with the virtual mass squared Q 2 of the partons in the shower, the infrared 
cut-off is denoted by Qq, and the splitting parameter z denotes the energy splitting. In this 
setup, the transverse momentum of the splitting is given by k\ ~ z(l — z)Q 2 for time-like 
and k\ ~ (1 — z)Q 2 for space-like splittings. The current implementation is limited to 
massless quarks, the extension to massive quarks is straightforward. 

Contrary to the way in which the Pythia 6 shower is implemented, the kinematics of 
each splitting is already fully constructed when the splitting is generated. This requires 
determining the daughter virtualities together with the splitting of the mother. For parton 
shower evolution in the vacuum this procedure is more involved than the traditional ap- 
proach of later correcting the kinematics of the splitting of the mother, but is advantageous 
for jet evolution in a medium. The reason for this is that in the medium one evolves not 
only in the shower ordering variable but also in time and going back to the mother splitting 
means going back in time, which is very cumbersome. 

Since the centre-of-mass energies of scatterings in the medium are typically rather 
small, it is unlikely that more than one initial state splitting is initiated in secondary 
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scatters. In this approximation 5 the partonic PDFs can be integrated analytically yielding 

Q 2 

S q (Q 2 , Ql)6(l - x) + J** S q (Q V) " s((1 -* )q2) P qq (x) , (3.3) 
Ql 

(3.4) 

Jfs^^ ^-^ P^), (3.5) 

Ql 

Q 2 

/|(x,Q 2 ) = |^S q (QV) as((1 P qg (x) , (3.6) 
Ql 

S g (Q\Ql)5{l-x), 

Ql 

3.4 Implementation of the LPM effect in JEWEL 

In Jewel it is assumed that the physical reasoning behind the probabilistic formulation 
of the LPM-effect is valid also outside the eikonal limit. Although the kinematics is very 
different, the procedure of adding scattering processes within the gluon formation time 
coherently by adding their momentum transfers and reweighting the actual emissions pro- 
ceeds in exactly the same way. 

A slight complication arises due to the fact, that - in contrast to the eikonal case - 
the formation time is not directly related to the momentum transfer but to the scale of the 
emission, which can be anything between the scale fW set by the momentum transfer and 
t c . The procedure for adding a momentum transfer coherently and determining the scale of 

(s) 

the associated radiation is outlined below. Here, t\ denotes the scale of scattering before 

fs) 

that last momentum transfer is added and t 2 is the scale including the last momentum 

fs) 

transfer. As the momentum transfers are added vectorially, the resulting t\ may be larger 

(s) . (s) 

or smaller than t\ . Before the last scattering, t\ was the starting scale of a parton shower 
and (provided > t c ) the parton shower may or may not foresee a gluon radiation at 
a scale t\ < Updating the starting scale set by the effective (coherent) scatterings to 
irj may increase or decrease the phase space for radiation and the existing emission (if 
there is one) has to be corrected for the corresponding change in phase space and emission 
probability. As a result an emission may have to be rejected, it may be assigned a new 
scale t2, or a new emission may be enforced where there was none before. All the different 
cases are summarised below. 

5 This approximation used here can easily be improved numerically, adding more emissions, if necessary. 



/ q q (x,Q 2 ) = 
f q i (x,Q 2 ) = 
/ g q (z,Q 2 ) = 

/|(x,g 2 ) = 

fi(x,Q 2 ) = 
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• 4 S) < t c : 

In this situation no radiation is possible, if there is an existing emission it has to be 
rejected. 

• existing radiation, t^' > t± . 

The probability for an emission is now larger and so the existing radiation cannot 
be rejected. The scale of the radiation, however, may have to be redetermined. 
The chosen t\ for the emission is kept with probability Sit^, t^), i.e. if there is no 
radiation between 4 and to it has to be in the old interval [i c ,ti S ^l. In the other 
case where we have radiation between t\ and t 2 a new scale t<i in this interval is 
chosen and the old one has to be discarded. 

(s) (s) 

• existing radiation, t\ > t\ > t c : 

In this case the probability for radiation has decreased and the existing emission has 
to be rejected with probability [1 - S(t { *\ t c )]/[l - S(t^\ Q]. If the emission is kept 
but its scale t\ is larger than t 2 a new scale t% in the new interval has to be chosen. 

• no existing radiation, > t^: 

The phase space for emission has increased and with probability 1 — S(t 2 \ 4 ) a new 

(s) (s) 

emission with scale ti to be chosen between t\ and t 2 is added. 

• no existing radiation, t± > t 2 > tc'- 

In this case nothing happens since the probability for radiation has decreased and no 
emission needs to be reweighted. 

The reweighting of the emission with 1/N sca t is straightforward and, as discussed in 
section 2.2, it will only be accepted if its final formation time is shorter than that of 
potential other emissions. 

3.5 Recoils, colour flows, and hadronisation 

In Jewel only the QCD evolution of jets and their interactions with the medium are 
simulated and not the complete event. The recoiling scattering centres can be traced and in 
principle they could be subjected to further interactions. This option is, however, disabled 
by default as these interactions are typically much softer than the jet-medium interactions 
and perturbation theory is not necessarily the correct language to describe them (and due 
to the rapid increase in number of interactions the run-time increases dramatically). As a 
way of estimating the effect of the jet on the medium and jet-background correlations the 
recoils may be kept in the event, removed or replaced with uncorrelated scattering centres. 

As in the vacuum parton shower every radiated gluon is a colour neighbour of the 
emitter. When the recoils are not traced they also don't exchange colour with the jet. 
When the recoils are kept in the event, they are assumed to be neighbours of their scat- 
tering partner and are inserted in the existing colour string. This obviously constitutes 
a simplification but avoids another complication. After all splittings and scatterings the 
event is hadronised using Pythia's implementation of the Lund string fragmentation. For 
the Lund string model to work properly it is essential that the strings are of sufficient 
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invariant mass. Naively trying to track all colour indices of all scattering centres without 
invoking any colour reconnection model indeed leads to large numbers of strings with very 
low invariant mass and the breakdown of the hadronisation model. 



3.6 Medium model 

Jewel can in principle be interfaced with any medium model that specifies the density and 
momentum distribution of scattering centres for any point in space-time. In order to have 
a simple, highly predictive model capturing all essential features of the medium without 
introducing too much additional, potentially intractable kinematics and, at the same time, 
allowing full control over the medium, we used a variant of the Bjorken model [65, 66] for 
this study. It describes the boost-invariant longitudinal expansion of an ideal quark-gluon- 
gas with three parameters only, namely the initial proper time T\ at which the hydrodynamic 
evolution sets in, the initial temperature T\ = T(t\) and the critical temperature T c of 
the deconfinement phase transition. The transverse profile is chosen such that the energy 
density is proportional to the density of wounded nucleons. The latter is, as all geometrical 
aspects, calculated in the framework of a simple Glauber model [67]. In practice this means 
that the initial temperature is translated into an initial energy density e; oc Tj 4 . The actual 
energy density profile at t; is then given by 

e(x,y,b,rO = e i ^f ( : y ' b) . with <n part > (b = 0) » -M- , (3.8) 

(ftpart)(o = 0) nR A 

where b is the impact parameter, n par t(x, y, b) is the density of participating nucleons in 
the transverse plane and we have for simplicity assumed a symmetric A + A collision. As a 
consequence the temperature is higher in the centre of the overlap region and it decreases 
with centrality. The (proper) time dependence is given by 

e{x,y,b,r) = e(x,y,b,n) ( — J and T(x, y, b, t) oc e 1/4 (x, y, b, r ; ) ( - J . (3.9) 

The time dependence of the particle density is then given by n(x,y,b,r) oc T 3 (x,y,b,r). 
We fix the critical temperature at T c = 165 MeV consistent with lattice results and since 
we are only considering interactions in the deconfined phase there are no scatterings when 
the local temperature has dropped below T c . 



4 Results 

The analysis of Monte Carlo events and all plots shown here were produced with Rivet [68] 
and the corresponding analysis codes contained within. For jets, typically the anti-/cj_ 
algorithm of [23] has been used throughout. 

4.1 Validation 

In the absence of a medium Jewel reduces to a standard (vacuum) parton shower, which 
was validated extensively against data from LEP and p + p collisions at LHC. 
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Figure 1. Comparison of Jewel+Pythia results to Lep data. LHS: Thrust distribution mea- 
sured by Aleph [69]. RHS: charged particle fragmentation function as a function of the scaled 
momentum x v = 2p\ l /^fs measured by Delphi [70]. 



In all results shown here, the strong coupling is running at one loop with a s (mz) = 
0.128, consistent with findings of other leading order parton shower calculations. The 
infrared cut-off of the parton shower was chosen as Qq = 1.5 GeV. As the Jewel par- 
ton shower is very similar to the Pythia virtuality-ordered shower and given the very 
reasonable description of the LEP data, a retuning of the hadronisation parameters was 
considered not to be necessary. 

As examples illustrating the generally very satisfactory performance of the parton 
shower (and hadronisation) at LEP the thrust distribution and the charged particle frag- 
mentation function are shown in Fig. 1. 

At the LHC the comparison of Jewel+Pythia results to the jet data is complicated 
by the underlying event, which is not included in the MC. This can, for instance, be 
seen in the inclusive (track) jet cross section shown in figure 2: While the agreement 
between the MC results and the data is satisfactory at high p±, the MC falls below the 
data at relatively small transverse momenta, where the contribution from the underlying 
event is largest. The same effect is visible in the (calorimetric) jet shape measurements 
(figure 3), where the Jewel+Pythia jet are significantly more collimated at low p± and 
again the agreement improves with increasing p±. A complementary observable, namely 
the azimuthal decorrelation shown in figure 4, is also described reasonably well. 

Figure 5 shows the comparison of Jewel+Pythia to the Atlas measurement of the 
charged particle fragmentation function in track jets. Here, the Jewel+Pythia frag- 
mentation consistently is too soft. The interpretation of this observation is, however, not 
straightforward, as the measurement is carried out for relatively low jet-p± and it is thus 
unclear to what extent there is a contamination from the underlying event. Two points 
hint at an issue in the MC in addition to possible underlying event contributions: Firstly, 
the charged particle density inside track jets shown in figure 6 shows better agreement 
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Charged jet cross section vs. p± (anti-/ct, R = 0.4) 
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Figure 2. Jewel+Pythia results for the inclusive jet cross section in different rapidity bins 
compared to Atlas data [71] in p+p at y/s = 7TeV. Jets are reconstructed using the anti-/cj_ 
algorithm with R=0.4 on tracks. 
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Figure 3. Jewel+Pythia compared to Atlas measurements of the jet shape p(r) [72] in p+p 
at y/s = 7 TeV. 



between MC and data, although also here the components missing at large r could be due 
to underlying event. Secondly, the jet shapes of calorimeter jets tend to be slightly too 




Figure 4. Jewel+Pythia compared to Atlas jet measurements of the azimuthal decorrela- 
tion [73] in p+p at */s = 7TeV. 




Figure 5. Jewel+Pythia results for the fragmentation function of charged particles in track jets 
for different jet-p_i_ bins compared to Atlas data [71]. 



collimated in the MC even at high p± (figure 3), which would be consistent with a softer 
fragmentation function. It should be kept in mind that the jet shape p[r) is an energy 
density, while p c h{f) is a particle density. 

The comparison of the Jewel shower to the measurement of the ratio of the 3-jet 
to the 2-jet cross section by CMS shown in Fig. 7 demonstrates that the parton shower 
provides a very reasonable description of the 3-jet matrix elements even well outside the 
collinear region. This confirms that leading order matrix elements in combination with a 
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Charged jet p c h(r) (anti-/cf, R = 0.4, y 0.0-1.9) 




Figure 6. Jewel+Pythia results for the charged particle density p c h{f) in track jets for different 
jet-p_i_ bins compared to Atlas data [71]. 
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Figure 7. The ratio R32 of the inclusive 3-jct to the 2-jet cross section as a function of the 
scalar sum of the transverse momenta Ht, the anti-fcj_ jets in this analysis have p± > 50GeV and 
R = 0.5 [74]. 



parton shower are a good approximation to radiative processes and thus provide a good 
estimate of the radiative energy loss. 

Fig. 8 shows a comparison with the neutral pion spectrum measured by Phenix and the 
charged particle spectrum measured by CMS, which form the baseline for the respective 
measurements of the nuclear modification factor Raa- In the former case, above p± ~ 



- 18 - 




Figure 8. LHS: Neutral pion spectrum in p+p collisions at y/s = 200 GeV simulated with 
Jewel+Pythia and compared to Phenix data [75]. RHS: Charged hadron spectrum in p+p 
collisions at y/s = 2.76 TeV simulated with Jewel+Pythia and compared to Cms data [18]. 



4 GeV, the Jewel+Pythia results agree with the data on a level of roughly 15 % over 
about 6 decades and thus provide a realistic baseline for the Raa determination. At LHC 
energies the hadron spectrum seems to be slightly softer in the simulation than in data. 
However, deviations of this size have only a very small effect on the nuclear modification 
factor. The uncertainty arising from this is much smaller than, for instance, the uncertainty 
arising from the choice of infrared regulator /iD- 

4.2 Hadron suppression 

To obtain a fair agreement with the measured nuclear modification factor at RHIC, 
Jewel+Pythia requires an initial temperature of Tj = 350 MeV at initial proper time 
n = 0.8 fm, see Fig. 9. This is remarkably consistent with the input parameters in fluid 
dynamic simulations of heavy ion collisions [77-79]. We note that at high p±, where the 
Monte Carlo results are reliable, they reproduce both the factor ~ 5 suppression and the 
approximately flat pj_-dependence seen in data. Varying the Debye mass by ±10% has a 
significant effect on the overall suppression, indicated by the grey band, but hardly affects 
the shape. The dependence on the exact choice of the regulator is clearly sizeable. The 
pragmatic approach is to regard it as a parameter, fix it together with the medium param- 
eters at one point (in this case the neutral pion suppression at RHIC) and use the same 
value for all calculations. But while this may be a practical working solution it does not 
address the underlying problem. 

The charged hadron multiplicities measured in heavy ion collisions constrain the initial 
entropy density of the system, sjTi cx dN/dy, s; oc e\/T\ oc Tj 3 and therefore allow to relate 
the initial temperatures at RHIC and at the LHC, 

nLHC _ t^rhic / 'i 

[IC dN/dy\ 

LHC 



1 i ~ J i \ _LHC jsr/x.l • \^- L 



T; 



dN/dy\ 



RHIC 
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Figure 9. Nuclear modification factor for neutral pions in Au+Au collisions at y/s = 200 A GeV 
in the 0-10% centrality class simulated with Jewel+Pythia and compared to Phenix data [76] 
(only systematic errors shown, statistical errors are smaller than the systematic ones everywhere). 
The grey band indicates a variation of the Debye mass by ±10 %. 

The observation of a factor 2.2 increase in the event multiplicity from RHIC to LHC is 
therefore consistent with an initial temperature T\ = 530 MeV at t\ = 0.5 fm at the LHC. 
There is some freedom in initializing the fluid dynamic evolution at the LHC at a different 
initial time T[, but this is numerically unimportant. At early times the parton shower is 
dominated by emissions at rather high scales initiated by the initial hard scattering, and 
this high virtuality protects the partons from medium-induced emissions and makes them 
insensitive to the medium at early times. Thus, the medium at the LHC is specified in terms 
of parameters fixed in Fig. 9. As seen from Fig. 10, the calculation of Jewel+Pythia 
then leads to a very good agreement with preliminary data of the nuclear modification 
factor at the LHC without any additional adjustments. This is a remarkable success of 
our model. 

To understand the characteristically different ^-dependencies of Raa at RHIC and at 
the LHC, we have performed control simulations in which a LHC-like distribution of hard 
processes is fragmented in a RHIC-like medium, and vice versa (data not shown). This 
showed that the significant change in the slope of -Raa(rl) from RHIC to the LHC can 
be attributed fully to the -y/s-dependence of the distribution of initial hard processes. We 
therefore conclude that a purely perturbative dynamics of parton energy loss supplemented 
by an arguably simple model of the medium whose characterisation matches physical ex- 
pectations, can account for the main features of the measured nuclear modification factors, 
including the strength of the suppression pattern, and its yfs- and pj_-dependence. 

At very large p± the nuclear modification factor continues to rise above unity. This is a 
purely kinematical effect that becomes visible at very large p± where the energy loss starts 
to vanish. The elastic scattering of energetic partons converts longitudinal into transverse 
momentum turning the pj_-spectrum harder. While this is a generic effect its size and 
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Figure 10. Nuclear modification factor for charged hadrons in Pb+Pb collisions at y/s = 
2.76 A TeV in the 0-5% centrality class simulated with Jewel+Pythia and compared to pre- 
liminary CMS [18] and Alice [16] data. 

turn-on point will to some degree depend on the medium model, as they are sensitive to 
the amount of scattering centres encountered in the forward direction. This effect would 
not be accessible in the standard eikonal treatment of parton energy loss that neglects the 
change in propagation direction of the projectile parton. We believe it is thus a generic 
testable prediction of our model. 

4.3 Modification of jets 

Analogously to the single-inclusive hadron suppression discussed in the previous section 
also the single-inclusive jet suppression has been measured by Alice and Atlas. In 
figure 11 and 12 the experimental results for different values of the jet radius are compared 
to Jewel+Pythia. The Alice results (which are relative to a pp baseline) and the Atlas 
results (which are relative to peripheral Pb+Pb spectra) for small jet radii are very well 
reproduced, but there is some discrepancy in the Atlas results for large jet radii. Also, 
in the R = 0.2 measurement by Alice, the pj_-dependence seems to be weaker in the MC 
than in the data. There is no indication of a significant rise of the jet suppression at large 
p± in the simulation, which shows that jets can behave differently from leading hadrons. 
The fact that the disagreement between Jewel+Pythia and the Atlas data increases 
with the jet radius R could hint at an issue with the background subtraction. 

For the simulation of the single-inclusive hadron spectra the recoiling scattering centres 
are kept in the event and hadronised together with the jets. In the case of jet analyses, 
however, the recoils have to be taken out, since the experimental jet reconstruction involves 
the subtraction of background. These two procedures are not exactly identical for several 
reasons. Firstly, in Jewel+Pythia the hadronic final state is not the incoherent sum 
of a jet component and a contribution from the recoils, due to the colour connections 
between the jet partons and the recoiling scattering centres. Secondly, in the experimental 
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Figure 11. Jewel+Pythia results for Raa of jets in Pb+Pb collisions at -^/snn — 2.76 TeV 
compared to Alice data [21] for two values of the jet radius (correlated systematic errors not 
shown) . 
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Figure 12. Jewel+Pythia results for Rcp of jets in Pb+Pb collisions at ^/snn — 2.76 TeV 
compared to Atlas data [22] for different values of the jet radius. The ratio is taken between the 
0-10% and the 60-80% centrality class. 



jet reconstruction procedure, only a background estimate based on the average activity in 
the event can be subtracted. While to some degree it is possible to deal with background 
fluctuations using unfolding techniques [80-82], the correlations between the jet and the 
background cannot be assessed in this way. There will thus always remain a residual 
uncertainty when comparing any MC results to jet data. Clearly, along the same line of 
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dijet asymmetry in Pb+Pb (^/s = 2.76 ATeV) 



azimuthal decorrelation in Pb+Pb (y^ — 2.76 ATeV) 




Figure 13. Di-jet asymmetry Aj = (-E_l,i — E±^)/(E±,i + E±,2) and azimuthal decorrelation in 
central Pb+Pb and p+p collisions at */Inn = 2.76 TeV, the cuts are the same as in the Atlas 
analysis [19]. 




Figure 14. Jewel+Pythia results for the fragmentation function D(z) in peripheral and central 
Pb+Pb events compared to Atlas data [26] (data points read off the plots, no errors shown). 



reasoning, the background subtraction becomes more involved as the jet radius increases 
and the observation that the discrepancy between Jewel+Pythia and data increases with 
jet radius may hint at a problem in this direction. 

Both Atlas [19] and CMS [25] have measured the di-jet pj_-asymmetry Aj and the 
azimuthal decorrelation in central Pb+Pb events. The Jewel+Pythia results obtained 
with the Atlas cuts are shown in Fig. 13. A direct quantitative comparison to either At- 
las or Cms data is not possible, as the data are not unfolded for the jet energy resolution, 
which is known to have a sizeable effect especially on the shape of the di-jet asymmetry. 
Qualitatively, however, the behaviour of the Jewel+Pythia results resembles observa- 
tions made in the data: There is a significant broadening of the Aj distribution in central 
Pb+Pb events as compared to the p+p baseline. The azimuthal decorrelation, on the other 
hand, hardly changes; there is only a mild increase at separations A<p ~ ir/2. 

The charged particle fragmentation function in jets in central Pb+Pb collisions as a 
function of the momentum fraction z and relative p± in central and peripheral events are 
shown in figure 14 and figure 15, respectively. Except for the region of very small z or p±, 
which is expected to be particularly sensitive to background modelling and subtraction, 
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Figure 15. Jewel+Pythia results for the fragmentation function D(p±) in peripheral and central 
Pb+Pb events compared to Atlas data [26] (data points read off the plots, no errors shown). 
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Figure 16. Jewel+Pythia results for the ratios of the fragmentation functions D{z) between 
central and peripheral Pb+Pb events compared to Atlas data [26] (data points read off the plots, 
only maximum of statistical and systematic errors shown). 



the agreement between the Jewel+Pythia results and the data is very reasonable given 
that for this measurement the background subtraction is more involved than for jets. It 
remains unclear, however, whether the observation, that the fragmentation function in p+p 
at yfs = 7TeV seems to be too soft, indicates that in Pb+Pb it is actually too hard. The 
ratio of the fragmentation functions in central and peripheral events shown in figure 16 
and figure 17 are less sensitive to such effects. Again, the agreement between data and 
MC improves for smaller values of the jet radius. In Jewel+Pythia the fragmentation 
function tends to become harder in more central events. This can easily be understood 
since the fragmentation of the hard core remains unaltered by the presence of the medium 
while the jet energy gets degraded. 
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Figure 17. Jewel+Pythia results for the ratios of the fragmentation functions D(p±) between 
central and peripheral Pb+Pb events compared to Atlas data [26] (data points read off the plots, 
only maximum of statistical and systematic errors shown). 



4.4 Uncertainties 

Uncertainties related to the basic assumptions and other aspects of the calculation (e.g. 
pdf uncertainties) are to some extent unavoidable, but can be quantified in our model. 
While a complete assessment of all possible uncertainties certainly is beyond the scope of 
this publication, we identify and discuss the most important sources of uncertainties and 
their respective sizes. 

Apart from the model of the medium the framework is essentially free of parameters. 
There are choices to be made concerning the infrared regulator of the matrix elements and 
the formation time, but once they are fixed they should not be changed. The infrared 
cut-off of the parton shower and the parameters of the hadronisation model can be tuned 
to LEP data and are not changed when going to nucleus-nucleus collisions. There is thus 
only very little room for adjustments and tuning. 

The effect of varying the infrared regulator /io, which presents the largest uncertainty, 
was already discussed in the context of the single-inclusive hadron suppression. In this 
section we show the effect of varying the formation time up and down by a factor 2 and 
using different PDFs (Cteq6l1[62]+Eps09lo[64], which is the default, Cteq6l1 without 
nuclear modifiactions and Mstw08Lo[83]+Eps09lo). Figure 18 shows the effect on the 
inclusive jet spectrum and the fragmentation function and figure 19 the di-jet asymmetry 
and azimuthal decorrelation. Varying the formation times leads to a ~ 20 % variation of the 
jet rate in central Pb+Pb collisions and a smaller but visible change in the di-jet asymmetry 
and azimuthal decorrelation while the fragmentation function is fairly insensitive. The 
uncertainty related to the choice of PDF is smaller and of the same size as the current 
statistical uncertainty on the MC results. 

In a similar way practically all aspects of the simulation can be varied and the uncer- 
tainties estimated. We find that the largest sources of uncertainties are the choice of the 
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Figure 18. Uncertainties due to varying the formation time (yellow band) and the pdf (blue band) 
on the inclusive jet spectrum (LHS) and the fragmentation function (RHS) in central Pb+Pb 
collisions. The statistical uncertainty on the default set-up is shown as error bars in the upper 
panel and as the red band in the ratio plots. 
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Figure 19. Uncertainties due to varying the formation time (yellow band) and the pdf (blue band) 
on the di-jet asymmetry (LHS) and the azimuthal decorrelation (RHS) in central Pb+Pb collisions. 
The statistical uncertainty on the default set-up is shown as error bars in the upper panel and as 
the red band in the ratio plots. 



infrared regulator and to a lesser degree the formation time. 
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5 Conclusions 



In this paper we report on the development of a novel description of jet quenching and 
its implementation into the Monte Carlo generator Jewel. The paradigm underlying the 
construction of our model is the language given by perturbative QCD, with the aim to 
enable a quantitative discussion of effects beyond it. In this context, hard interactions 
of a jet with the medium resolve its constituents as quasi-free partons and should be de- 
scribed in perturbation theory 6 . In reinterpreting induced radiation as regular, although 
infra-red-regulated, 2 — > 2 parton scatterings, supplemented with a parton shower we ar- 
rive at a framework whose dynamics is entirely based on standard perturbative technology 
also used in the simulation of p + p collisions and which is minimal in its assumptions. 
It is formulated entirely in general kinematics 7 and thus overcomes limitations of previ- 
ous approaches which were entirely based on results obtained in eikonal or close-to-eikonal 
kinematics [48]. Consequently, there is no distinction between elastic and inelastic scat- 
tering. Rather, both are possible outcomes of the same process and the parton shower 
cut-off defines the separation between the two. Also the distinction between vacuum and 
medium-induced radiation has become obsolete, as they are fundamentally the same. The 
interplay between all sources of radiation is governed by the formation times and generally 
a certain emission cannot be classified as vacuum-like or medium-induced in a meaningful 
way. 

The emerging framework is well constrained by data from elementary reactions and 
does not leave much room for tuning. Apart from the model of the medium there is 
only a single quantity with a significant uncertainty that should perhaps be regarded as 
a parameter, namely the infrared regulator of the matrix elements. All uncertainties re- 
lated to assumptions or other aspects of the calculations can be quantified, examples for 
the numerically most important uncertainties have been presented here, but a complete 
characterisation of all possible uncertainties is beyond the scope of this publication. 

Our new framework does not make any assumptions about the medium other than that 
is partonic. The phase space distribution does not matter and in principle any distribution 
can be interfaced. The results presented here were obtained with a simple Bjorken type 
model. Nevertheless, after fixing the parameters with RHIC data, the agreement with very 
different jet and leading hadron measurements at the LHC is very reassuring. This leads 
us to conclude that our calculation, although certain aspects are simplified, captures the 
relevant physics underlying jet quenching. At the same time the overall good description of 
data lets us speculate in how far more complicated medium models may obfuscate simple 
physical realities. This clearly is a question that could only be answered seriously by 
interfacing more medium models and cataloguing the results obtained with them. 

However, we would like to stress that for some observables there appear to be unre- 

6 In general, the scale above which the components of a final state parton shower can resolve partonic 
constituents of the medium as quasi-free partons depends on properties of the medium [84]. Comparing a 
perturbative formulation of jet quenching to data can thus give access to the nature of medium constituents. 

7 Although the present code goes far beyond the eikonal kinematics, it is clear how to recover this limit 
of the full emission pattern encoded here. 
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solved and potentially unresol valve issues with the subtraction of backgrounds, where the 
procedure in the analysis of MC events differs from the experimental one and may lead to 
discrepancies. The observation that generally the agreement between MC and data gets 
worse for larger jet radii points in this direction. Also, it should be noted that there is a 
class of measurements that are sensitive to the back- reaction of the jet on the medium and 
that therefore our calculation cannot be expected to describe. This problem could only be 
overcome by trying to extend our model to the fully inclusive treatment of the full collision 
- something clearly beyond the scope of this paper. 

Apart from the obvious option of testing different models for the medium we have 
hinted at various other ways of how our calculation can systematically be improved using 
a well-understood perturbative language. 
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